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Spirals and clumps in V960 Mon: signs of planet formation via gravitational instability around an FU Ori star? 
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ABSTRACT 


The formation of giant planets has traditionally been divided into two pathways: core accretion and gravi- 
tational instability. However, in recent years, gravitational instability has become less favored, primarily due 
to the scarcity of observations of fragmented protoplanetary disks around young stars and low occurrence rate 
of massive planets on very wide orbits. In this study, we present a SPHERE/IRDIS polarized light observation 
of the young outbursting object V960 Mon. The image reveals a vast structure of intricately shaped scattered 
light with several spiral arms. This finding motivated a re-analysis of archival ALMA 1.3 mm data acquired 
just two years after the onset of the outburst of V960 Mon. In these data, we discover several clumps of con- 
tinuum emission aligned along a spiral arm that coincides with the scattered light structure. We interpret the 
localized emission as fragments formed from a spiral arm under gravitational collapse. Estimating the mass of 
solids within these clumps to be of several Earth masses, we suggest this observation to be the first evidence of 
gravitational instability occurring on planetary scales. This study discusses the significance of this finding for 
planet formation and its potential connection with the outbursting state of V960 Mon. 


Keywords: FU Orionis stars (553) — Gravitational instability (668) — Observational astronomy (1145) — Planet 
Formation (1241) 


1. INTRODUCTION The core accretion scenario, characterized by the continu- 
ous growth of dust particles followed by runaway gas accre- 
tion, is the prevailing formation scenario for gas giant plan- 

philipppweber@ gmail.com ets (Pollack et al. 1996). However, when it comes to directly 
imaged giant planets and brown dwarfs found at significant 
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distances from their host stars, most models face difficulties 
in producing sufficiently massive objects within the expected 
lifetimes of gaseous disks (e.g. Emsenhuber et al. 2021). As 
a result, gravitational instability (GI) emerges as a prominent 
alternative for planet formation in these regions (see Kratter 
& Lodato 2016, for a review), as it is believed to operate in 
the outer protoplanetary disk. 

The underlying physical mechanisms that drive the forma- 
tion of massive gas giant planets through gravitational insta- 
bility in circumstellar disks may also be responsible for the 
Observed episodic accretion events in young stellar objects 
(Armitage et al. 2001; Fischer et al. 2022). According to 
this model, protostars undergo intermittent and transient, yet 
highly efficient accretion phases, known as “FUor events” 
(named after the prototype event observed in FU Orionis). 
Those events are characterized by a significant increase in ac- 
cretion luminosity occurring over an annual timescale. FUor 
objects pose an excellent laboratory to study planet forma- 
tion, as the sudden increase in brightness of the central source 
illuminates and heats the surrounding environment. This has 
been exploited in the past for the first direct detection and 
characterization of the water ice line (Cieza et al. 2016; To- 
bin et al. 2023). For an extensive discussion of FUor objects, 
existing observational evidence, and its theoretical interpre- 
tations we refer to dedicated review articles (Audard et al. 
2014; Fischer et al. 2022). 

V960 Mon (2MASS J06593158-0405277) is a bona fide 
FUor object (Connelley & Reipurth 2018; Cruz-Sáenz de 
Miera et al. 2023) that has been in an outbursting state since 
2014 (Maehara et al. 2014). The distance to the object has re- 
mained a subject of controversy over the years. Initially, Kim 
et al. (2004) proposed a kinematic distance of 2.3 kpc based 
on CO observations and association with the molecular cloud 
S 287. However, when compared to stellar evolution models, 
a distance of approximately 450 pc was estimated (Kóspál 
et al. 2015). More recently, Kóspál et al. (2021) used a dis- 
tance of 1574 pc from the Gaia DR2 dataset. These contrast- 
ing results underscore the sensitivity of the inferred distance 
to the chosen method of measurement. 

Kóspál et al. (2015) analyzed several pre-burst archival 
datasets across the spectrum and proposed eight compan- 
ions in the local field to be pre-main sequence objects in a 
T Tauri stage. From this, the authors suggest that V960 Mon 
is not isolated. Most recently, Cruz-Sáenz de Miera et al. 
(2023) observed the blue- and red-shifted parts of an out- 
flow in !?CO emission approximately along the line-of-sight. 
From !?CO measurements, they estimate a massive envelope 
of ~ 0.6 Mo. 

In this Letter, we present and analyze prominent large- 
scale spirals around the FUor V960 Mon observed with 
SPHERE/IRDIS in polarized light. By using archival ALMA 
1.3 mm continuum data we reveal dust clumps within these 


structures. We describe the observational setup in Section 2, 
present the results in Section 3, and discuss their implica- 
tions in Section 4. These observations potentially connect GI 
clumps to a recent FUor outburst and offer an unprecedented 
opportunity to characterize gravitational instability at planet- 
forming scales. 


2. OBSERVATION AND DATA REDUCTION 


We observed V960 Mon in H-band (A = 1.625 um) in the 
night of 2016 December 17 (programme-ID: 098.C-0422(B), 
PI: L. Cieza) in dual-beam polarimetric imaging mode (DPI, 
de Boer et al. 2020; van Holstein et al. 2020) with the In- 
fraRed Dual-band Imager and Spectrograph (IRDIS, Dohlen 
et al. 2008) of VLT/SPHERE (Beuzit et al. 2019). The tar- 
get was part of a SPHERE survey to detect the scattered 
light from the environments around six episodic accreting 
FUors and EXors, which will be presented in Zurlo et al. 
(in prep). The observing conditions during operation were 
excellent (seeing ~ 0/4, vwina 73.5 ms-!, ro ~10 ms). We 
used the N. ALC. YJH.S coronagraph (185 mas diameter, Car- 
billet et al. 2011; Guerri et al. 2011) centered on V960 Mon 
and took a total of four half-wave plate (HWP) cycles of four 
object frames each. Each object frame was exposed to a de- 
tector integration time (DIT) of 64 s, adding up to a total inte- 
gration time of 17 min. Immediately before and after taking 
the polarimetric science frames, we took star center frames 
(DIT of 64 s) to accurately infer the stellar position behind 
the coronagraph and star flux frames (DIT of 2 s) to calibrate 
the observed intensities. 

We used the IRDAP reduction pipeline (IRDIS Data re- 
duction for Accurate Polarimetry, version 1.3.4, van Hol- 
stein et al. 2020) to process the SPHERE/IRDIS data and 
extract the Stokes Q and U components of linear polariza- 
tion. The pipeline incorporates a comprehensive model of 
the SPHERE optical system, enabling direct correction for 
instrumental polarization and polarization crosstalk without 
relying on the data. IRDAP automatically corrects for the re- 
duction of the unresolved stellar polarization due to the cen- 
tral star. We further subtract the unresolved polarization car- 
ried in the stellar halo of the close-in companion according to 
Weber et al. (2023, Appendix B). The total linearly polarized 
intensity is then calculated as: 


PI = YQ? +U?. (1) 


We refrain from using the conventional Q,/U, representa- 
tion commonly employed in polarized light imaging, as it 
is not suitable for capturing the polarized intensity if off- 
centered light sources are present within the local environ- 
ment. 


3. RESULTS 


Figure 1 shows the IRDIS polarized light H-band image 
for V960 Mon. The frame is centered on the primary star 
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Figure 1. SPHERE/IRDIS linearly polarised intensity (PJ) obser- 
vation of V960 Mon in H-band. PJ is normalized to its maximum 
value (6.6 mJy arcsec~*) and on a logarithmic scale. The frame is 
centered on the primary, which was covered by the coronagraph dis- 
played by a black circle. The projected scale is shown in the lower 
right assuming the distance of 2189 pc. 


which is concealed behind a coronagraph marked by a black 
circle. The image displays a vast S-shaped structure of scat- 
tered light extending along the north-south axis. Both the 
northern and southern parts are comprised of at least two 
adjacent spiral arms each. Assuming the Gaia DR3 dis- 
tance (d = 2189 + 281 pc, Gaia Collaboration et al. 2022), 
the projected extent of these spiral arms is several thousand 
au. Although we should approach the Gaia DR3 distance 
with caution due to the relatively large renormalized weight- 
ing unit error (RUWE) of 3.75, it is noteworthy that the 
Gaia DR3 distance calculated for the assumed companion 
UCACA 430-024261 (also referred to as V960 Mon N) is 
similar at 2606+346 pc (Gaia Collaboration et al. 2022), with 
a sufficiently small RUWE of 0.93, indicating a reliable as- 
trometric solution. However, it is important to interpret the 
projected spatial scales presented in this study while consid- 
ering the controversy associated with the distance measure- 
ment of the object (as introduced in Section 1). 

We confirm a close stellar companion southeast of the 
coronagraph which appears both in polarized and unpolar- 
ized intensity. This object was previously detected in J- and 
K-band, with Pa’ and Bry emission at a distance of 227 mas 
and a PA of 131.4 deg (Caratti o Garatti et al. 2015). Here, we 
measure the companion to be at a distance of 237--4 mas and 
at a PA of 136.7+1.0 deg with respect to the primary. No- 


tably, this companion is co-located with the scattering struc- 
tures, meaning there is no evidence of its orbit being cleared 
(see Fig. 1). The detection of the companion in both polar- 
ized and unpolarized intensity implies the presence of small 
dust grains in its immediate vicinity or indicates significant 
scattering along the line of sight. 

Motivated by the spectacular scattered light image seen in 
Fig. 1, we conducted a re-examination of archival ALMA 
band 6 data (programme-ID: 2016.1.00209.S, PI: Takami) 
previously published in Kóspál et al. (2021). We describe 
our data reduction and the differences to Kóspál et al. (2021) 
in Appendix A, with an exploration of the imaging parame- 
ters to test for the robustness of the observed structures. 

In the left panel of Figure 2 we show the ALMA 1.3mm 
continuum image and compare it to the polarized light image 
in the right panel. Both panels superimpose ALMA contours 
at 3, 4, and 5 Oyms levels (with o,44, = 28 uJy beam !). Re- 
markably, these contours reveal multiple spatially-separated 
fragments of continuum emission, reaching up to 7c. They 
align along a clockwise-opening spiral arm originating from 
the primary source. The left panel of Fig. 2 shows that these 
clumps roughly coincide with the southern spiral structure 
Observed in scattered light. We note that although the respec- 
tive observations were taken within a small temporal sepa- 
ration, a perfect alignment is not expected as the scattered 
light image probes the illuminated surface of the structures, 
while the ALMA continuum image traces emission from dust 
in cold, dense clumps that are likely optically thick in the 
near-infrared (NIR). In Appendix B, we show that the mm 
clumps can align with the scattered light spiral by assuming 
an adequate orientation for the spiral and a function for the 
scattering surface. 

Both the emission centered on V960 Mon and the individ- 
ual clumps remain unresolved by the beam of the observation 
(0714 x 0720 for natural weighting). We note that the close 
companion seen in the SPHERE data does not exhibit any 
significant counterpart in millimeter emission. 

Assuming the 1.3mm continuum emission is optically 
thin, we can estimate the clumps' masses based on a given 
temperature and dust opacity. Considering a typical opac- 
ity of K1.3mm = 2.3cm?g~! (Beckwith et al. 1990) and a 
temperature of 50 K, the clumps correspond to solid mate- 
rial masses ranging from 3 to 10 Mg, which corresponds to 
a gas mass of 1 to 3 Mjup, assuming a gas-to-dust mass ratio 
of 100. However, it is important to note that both the opac- 
ity and temperature can significantly deviate from these val- 
ues, depending on dust properties and local thermodynam- 
ics. Also, the gas-to-dust value can be significantly smaller 
where dust accumulates. The calculated masses exhibit an 
approximate inverse dependence on the assumed temperature 
and show a quadratic proportionality with the assumed dis- 
tance to V960 Mon. This indicates that if the previously as- 
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Figure 2. Comparison of ALMA band 6 continuum image in mJy beam~* using natural weights (left panel) with SPHERE/IRDIS polarized 
light image (right panel, see Fig. 1). Both images are overlaid with contours of the ALMA continuum at levels of 3, 4, and 5 ous (with orms = 
28uJy beam" +). The clean beam of the ALMA observation is shown in white in the bottom right corner of the left panel (0/14 x 0"20).The 


arrows in the left panel mark the dust clumps. 


sumed distances were used, the inferred clump masses would 
be even smaller. Consequently, we argue that the continuum 
emission effectively traces clumping occurring at scales rel- 
evant to planet formation. 


4. DISCUSSION 
4.1. Origin of spiral structures 


Spiral structures have been observed in several protostellar 
systems (see table 2 in Bae et al. 2022). These spirals exhibit 
variations in size, number, contrast, and pitch angle, suggest- 
ing different physical origins. Commonly discussed possi- 
bilities in the context of protoplanetary disks are a massive 
external companion (e.g. Dong et al. 2015) or vortex (van 
der Marel et al. 2016; Huang et al. 2019), gravitational insta- 
bility (e.g. Lodato & Rice 2005), an inner binary (Price et al. 
2018), a stellar fly-by (e.g. Clarke & Pringle 1993) or combi- 
nations of these processes (Thies et al. 2010; Pohl et al. 2015; 
Meru 2015). 

The complex environment surrounding V960 Mon 
presents a challenge for pinpointing the precise cause of the 
scattered light structure. The significant measured envelope 
mass (Cruz-Sáenz de Miera et al. 2023) implies a further 
potential association between the large-scale spirals and in- 
falling material (Lesur et al. 2015; Hennebelle et al. 2017; 
Kuffmeier et al. 2018). Another plausible explanation could 
be the capture of a close-by cloudlet (Dullemond et al. 2019). 


Additionally, the presence of multiple objects in the immedi- 
ate vicinity (Kóspál et al. 2021) introduces the possibility of 
considering them as potential candidates of a past fly-by. 
Kóspál et al. (2021) estimate that the mass correspond- 
ing to the unresolved emission around the primary could 
be as high as 0.33 Mc (updated from the original value of 
0.17 Mo considering the Gaia DR3 distance). The authors 
calculate that a disk of such a mass around an approximately 
solar-mass star would be susceptible to gravitational instabil- 
ity beyond a critical radius. However, simulations indicate 
that disks undergoing gravitational instability typically can- 
not maintain spiral arms beyond a radius of 100au for long 
periods of time, as the disk tends to fragment at larger radii 
(Rafikov 2005; Cossins et al. 2010; Zhu et al. 2012). 


4.2. Gravitational Fragmentation 


For V960 Mon, it remains uncertain which of the preced- 
ing formation scenarios accurately describes the environment 
surrounding the dust clumps. So far, Keplerian rotation was 
only detected in the optical and NIR (Park et al. 2020), where 
the molecular line profiles trace material much closer to the 
star. We speculate that the spirals are located between the in- 
ner envelope and outer disk and refer to disk equations with 
reservation. 

To our knowledge, the only process capable of explaining 
the fragmentation of a spiral arm into clumps (as witnessed 
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in Fig. 2) is the gravitational instability. The initial onset of 
gravitational instability is governed by the Toomre instability 
criterion, with Q « 1, where Q is defined as 


cs 
TGE’ 
with the local sound speed cs, the Keplerian frequency Q, the 
gravitational constant G and the surface density X. During 
gravitational instability, a disk is expected to generate large- 
scale spiral arms (Zhu et al. 2012) strikingly similar to those 
seen in the scattered light around V960 Mon in Fig. 1. Those 
spirals are expected to induce shock waves throughout the 
disk, which heat the disk material and regulate or potentially 
prevent further gravitational collapse. The ultimate fate of 
the gravitational collapse hinges on the efficiency of the disk 
material in radiating its thermal energy, directly relating to 
the cooling timescale (Gammie 2001; Rafikov 2005). 

The theoretical/numerical works of Takahashi et al. (2016) 
and Brucy & Hennebelle (2021) have predicted that if cool- 
ing is efficient enough, fragmentation of a spiral arm can oc- 
cur as the second stage of a two-step gravitational instability 
process. This is characterized by an adjusted instability cri- 
terion of Q < 0.6 within the spiral arm. 


Q= (2) 


4.3. Fate of clumps 


The detection of clumps within the ALMA band 6 observa- 
tion of V960 Mon marks a significant milestone by providing 
the first concrete evidence of a fragmenting spiral arm. This 
discovery indicates the involvement of gravitational instabil- 
ity in the formation of planetary-mass clumps and the evo- 
lutionary processes occurring within protoplanetary disks, at 
least in certain instances. 

What will happen to these clumps in the long run? Nu- 
merous studies have examined the fate of fragments resulting 
from gravitational instability. A comprehensive overview can 
be found in Kratter & Lodato (2016). In this context, several 
key concepts should be highlighted. 

One possibility is that the clumps disintegrate shortly af- 
ter their formation. This can occur due to tidal interactions 
with other clumps, rapid encounters with spiral structures, or 
insufficient cooling associated with tidal destruction. 

Furthermore, Vorobyov & Basu (2015) discovered that 
clump interactions can lead to rapid accretion onto the cen- 
tral star, potentially causing an accretion outburst. As a re- 
sult, a secondary clump may be ejected into the interstellar 
medium. These ejected clumps may be precursors of free- 
floating planets (Sumi et al. 2011) or brown dwarfs (Basu & 
Vorobyov 2012). 

Moreover, Zhu et al. (2012) conducted numerical simula- 
tions demonstrating that the fates of the clumps depend on 
the migration speed, cooling, and accretion efficiency. Slow 
migration leads to the formation of a massive companion 
opening a gap, while fast migration leads to tidal destruction. 


Nayakshin (2010) showed that in the tidal destruction sce- 
nario, the outer, more volatile components of the clumps are 
sequentially stripped away, and can ultimately leave a solid 
core in the inner region (a process labeled tidal downsizing). 

Finally, the interest in the formation of planetary cores 
through gravitational instability has been rekindled with the 
inclusion of solid material in the process of gravitational in- 
stability (Baehr et al. 2022). The authors found that un- 
der conditions specific to different dust sizes, overdensities 
can collapse and survive to give rise to planetary embryos. 
This holds the potential to significantly reduce the timeline 
of planet formation and offers a promising explanation for 
the detection of planetary signposts within the outer disks of 
even very young systems (Baehr 2023). 


4.4. FUor event 


Additional evidence supporting the presence of GI can be 
attributed to V960 Mon's classification as a bona fide FUor 
object. It has long been anticipated that FUor objects would 
serve as promising candidates for detecting indications of GI 
in their surrounding disks. Several studies have demonstrated 
that episodic accretion events can be triggered by the inter- 
play between GI and MRI (Armitage et al. 2001; Zhu et al. 
2009; Martin & Lubow 2011). Vorobyov & Basu (2005) fur- 
ther established that the inward spiraling clumps resulting 
from GI can lead to episodes of intense accretion, mirroring 
the observed behavior of FUor objects. 

The environments of most FUor objects look disrupted 
when observed the NIR scattered light (Liu et al. 2016, 
Takami et al. 2018, Zurlo et al. in prep). The mm- 
continuum, however, had thus far exhibited no signatures of 
any dynamical perturbances (see V883 Ori, Cieza et al. 2016; 
V900 Mon, Takami et al. 2019; FU Orionis, Pérez et al. 2020; 
HBC 494, Nogueira et al. 2023). This raises an immediate 
question regarding the persistence of the fragmented struc- 
tures for mm grains. According to Klahr & Schreiber (2020), 
the timescale for a dust clump's contraction is estimated to 
be tT ~ (95,Q)-!, where S, represents the Stokes num- 
ber representative of the dust's dynamical behavior. For dust 
particles of S, = 1 at 100au around solar-mass stars, this 
contraction timescale would be as short as ~ 20 years. After 
contraction, the clumps might not be detectable any longer. 


5. CONCLUSIONS & IMPLICATIONS 


We presented a polarized light image obtained with the 
VLT/SPHERE instrument, revealing the intricate system sur- 
rounding the FUor object V960 Mon. The image exhibits 
remarkable spiral arms extending over hundreds of astro- 
nomical units. Building upon these findings, we reanalyzed 
archival ALMA data, which provided deeper insights into the 
structure of the spirals, uncovering clumpy features located 
slightly offset with respect to the scattered light emission. We 
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quantified the properties of these clumps and estimated their 
dust and gas masses, marking the first detection of clumps 
in the planetary mass regime. The observed characteristics 
of V960 Mon closely resemble the predictions from simula- 
tions of gravitational instability, underscoring the suitability 
of FUor objects as laboratories for studying planet formation. 
The timing of the ALMA observations is particularly sig- 
nificant, as they were conducted merely two years after the 
stellar outburst. However, it is crucial to acknowledge that 
the FU Ori phenomenon might involve multiple triggering 
mechanisms, with gravitational instability-induced fragmen- 
tation representing just one possibility among several others. 
To gain further insights, a follow-up of FUor outbursts de- 
tected by the Legacy Survey of Space and Time (LSST) using 
the Vera C. Rubin telescope, specifically targeting millimeter 
clumps, could help establish the prevalence of such features 
around FUor sources in the early stages of outburst events. 
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Figure 3. Applying synthesis imaging reconstruction techniques to the ALMA band 6 data. a) and b) show the tclean reduction, a) using 
natural weighting (same as in Fig. 2) and b) using briggs-weighting with robust=0.5. The resulting clean beam sizes are 0714 x 0"20 and 
010 x 0716, respectively. Panels c) and d) show the GPUVMEM reconstruction technique applied to the data. Panel c) was restored using 
the clean beam of panel a), panel d) shows the image model. In all panels, the clean beam is displayed in translucent white in the bottom 
right corner. We measure the rms-noise, o for each image individually, written in the top left corner of each respective image and in units of 
wy beam ^! . The contours are shown at levels of 3c, 4c, and 5c. The position of the central star is marked by a +. The GPUVMEM model does 
not allow a measurement of the rms-noise. 


APPENDIX 


A. ALMA DATA REDUCTION 


Observations of V960 Mon were taken as part of the project 2016.1.00209.S (PI: M. Takami) and previously published in 
Kóspál et al. (2021). The data were obtained using a combination of ALMA configurations; these included an extended 12-m 
array configuration (2017 July 27, maximum baseline 3.7 km, pwv 0.47 mm) and a compact configuration (2017 April 20, 460 m, 
2.2 mm). Two spectral windows of bandwidth 1.875 GHz were dedicated to the continuum at central frequencies 216.877 GHz, 
and 232.178 GHz. Three other spectral windows were set up for molecular line observations, each with a total bandwidth of 
0.059 GHz, positioned with central frequencies 230.514 GHz, 220.375 GHz, and 219.537 GHz to cover the rotational 2-1 transi- 
tions of “CO, !?CO, and C180 respectively, using channel widths of 15.259 kHz, 30.158 kHz, and 30.158 kHz. We reduced the 
CO molecular line data to examine potential dynamical features evident in the moment 1 maps. However, the analysis yielded 
inconclusive results due to the limited on-target observation time of only — 12 min in total. Data calibration was carried out using 
the CASA pipeline in software version 5.6.1-8. Imaging and further analysis use CASA version 6. In our Letter, we present 
continuum imaging that combines all frequency channels devoid of line emission, employing a total aggregate bandwidth of 
3.445 GHz. 

We utilized the standalone version of the automated self-calibration module for the ALMA Pipeline (Tobin et al., in prep.)! 
to perform data self-calibration. The self-calibration process was performed separately for each array configuration, with two 
iterations of phase-only self-calibration. Following self-calibration, the compact configuration data showed an increase in signal- 
to-noise ratio of 17%. However, the more extended baseline dataset exhibited only a marginal improvement ranging from 1% to 
2%. Subsequently, both datasets were combined by concatenation to generate the final composite image. 

The imaging presented within this Letter in Fig. 2 applies natural weighting to the visibilities, resulting in a synthesized beam 
size of 0/14 x 0720. To provide a comparison, we also applied an alternate weighting scheme of briggs (robust=0.5) in Figure 3, 
which improved the beam size to 0710 x 0"16. However, this enhancement comes at a slight cost to the signal-to-noise ratio of 
the clumps. In Fig. 3, it is evident that the size of the clumps decreases as the beam size becomes smaller, indicative of their 
approximate point source nature. 

To ensure that the identified clumps are not artifacts resulting from the tclean image synthesis process, we reconstructed 
an independent image with the maximum entropy method (MEM) using the GPUVMEM? package (Cárcamo et al. 2018). To 


! https://github.com/jjtobin/auto_selfcal 


? https://github.com/miguelcarcamov/gpuvmem 
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Figure 4. Fitting consistent spirals to the images: The left panel shows the ALMA data with the 3, 4, and 5 orms contours. The white crosses 
represent the data points used for tracing the spiral arm. The black curve corresponds to a function fitted to the deprojected data points, which 
are then re-projected onto the image plane. The SPHERE observation in the right panel includes the same white crosses delineating the mm- 
spiral. Additionally, the yellow crosses represent the projection of the white crosses onto the scattering surface, assuming the function for the 
scattering height, hacac, described in the text. 


apply this method, we utilized a penalization factor of A = 0.01 in the entropy regularization term (see Cárcamo et al. 2018, 
for a description of the relevant parameters). The MEM reconstruction is an independent and unsupervised algorithm known for 
its ability to provide the highest resolution while maintaining sensitivity. Since the MEM algorithms forward-model visibility 
measurements using nonparametric models of the sky brightness distribution following well-informed image priors, the model 
image provided by GPUVMEM is a good representation of the actual sky brightness. In order to account for residuals and have an 
image comparable to the one produced by t c1ean, we also present a "restored" version of the model, which convolves the model 
with the natural weighted beam and adds the residuals. In Fig. 3, panel c) exhibits the restored image that matches the resolution 
of the natural weights. Remarkably, the restored version closely resembles the t c1ean image, underscoring the independence 
of the identified features from the image synthesis process. Fig. 3 panel d) displays the image model that effectively enhances 
resolution to approximately one-third of the original scale. The model image suggests that the clumps are significantly smaller 
than the natural beam size, corroborated by the contours in panel b). Both panel b) and d) illustrate that the emission centered on 
the primary source in panel a) encompasses two additional clumps that merge with the central feature for the larger beam sizes. 

A final word about the previous publication of the ALMA data: the clumps are not prominent in the same ALMA data presented 
in Kóspál et al. (2021) likely due to the authors' primary focus on resolving the inner disk and their use of uniform gridding in 
the imaging process. This approach compromises the signal-to-noise ratio, resulting in the faint clumps being closer to the noise 
level. However, it should be noted that the clumps also exhibit partial 5c-level signals when applying uniform weights (see fig. 1 
in Kóspál et al. 2021). In our study, we were guided by the IRDIS observation, motivating us to prioritize the detection of subtle 
features. Therefore, we utilized natural weighting in our imaging process to optimize for point-source sensitivity. This decision 
helped us highlight and investigate the previously unnoticed fragmenting spiral in greater detail. 


B. FITTING THE SPIRAL 


In Section 3, we discussed that the traced spiral arm of the mm-clumps in the ALMA image does not precisely match the spiral 
arm observed in the SPHERE image. We commented that this discrepancy arises due to the distinct nature of the observations. 
ALMA primarily captures thermal radiation emitted directly from the dust grains, while SPHERE images the surface where the 
dust structure becomes optically thick to incident stellar irradiation and scatters those photons to the observer. Consequently, the 
Observed radiation originates from two different planes, and when observed under an inclination, the projection of the spiral can 
lead to different locations in the image plane, resulting in a incongruent visual effect. 
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Our objective is to investigate whether and under what conditions the mm-clumps can be aligned with the closest spiral arm 
observed in scattered light. To achieve this, we make the assumptions that all mm-clumps lie within a single plane and that the 
scattered light arises directly above the mm-spiral. The comparison is carried out in two steps. First, we deproject the imaged 
ALMA spiral by assuming fixed values for inclination and PA. We fit an exponential spiral function of the form ¢(r) = a exp(br) 
to the deprojected ALMA data points, where r represents the deprojected radial distance to the center and ¢ denotes the azimuthal 
coordinate in the deprojected plane, increasing in the counterclockwise direction. The R?-value of the fit is best for an inclination 
of 63? and a PA of 120?. However, it is important to note that there were several other combinations of inclination and PA that 
resulted in nearly equally satisfactory outcomes. The spiral parameters at this deprojection yield an amplitude of a — 355? and 
a rate of change of b = 1.2 x 1074 au^ !. In the second step, we explore the parameter space of the scattering surface function, 
hscat = ho x (r/ ro)”, to identify the values that align the mm-spiral with the scattered light spiral. We set ry = 500 au and find 
that the best results are obtained for hy = 0.15 and v = 0.6. 


